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Correlations in the hadron distributions produced in relativistic Au+Au collisions are studied in 
the discrete wavelet expansion method. The analysis is performed in the space of pseudorapidity 
(I??! 1) a-nd azimuth (full 2ti) in bins of transverse momentum (pt) from 0.14 < pt < 2.1 GeV/c. In 
peripheral Au+Au collisions a correlation structure ascribed to minijet fragmentation is observed. 
It evolves with collision centrality and pt in a way not seen before which suggests strong dissipation 
of minijet fragmentation in the longitudinally-expanding medium. 

PACS numbers: 24.60.Ky, 25.75.Gz, 25.75.-q 



The study of the bulk properties of strongly interact- 
ing matter under extreme conditions at the Relativistic 
Heavy Ion Collider (RHIC) is producing a number of tan- 



talizing results 0. The physics of central Au+Au colli- 
sions at RHIC is clearly much more complex than a mere 
independent superposition of nucleon-nucleon collisions, 
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while the issues of possible collectivity and of the degree 
of "thermalization" of the bulk hadronic medium remain 
open. Substantial equilibration, especially in a short- 
lived finite system, may imply that during the evolution, 
there was a large number of degrees of freedom involved, 
such as would occur in a partonic medium or quark- 
gluon plasma Q. Equilibration in heavy ion collisions 
has been studied via its effects on single particle spectra 
and identified particle ratios. It progressively erases cor- 
relations, starting with the smallest features 0,0. Sur- 
viving correlations produced by hard scatterings early in 
the collision provide a sensitive monitor of the degree of 
equilibration of the medium. In contrast, traversal of the 
QCD phase boundary may create specific dynamical cor- 
relations . Therefore correlations observed in the final 
state are potentially affected by competing mechanisms. 
This makes the question of equilibration a quantitative 
one and warrants a study of correlations among the ma- 
jority of hadrons over a range of momentum scales. This 
Letter reports such a study. 

In high energy elementary collisions, hadrons originate 
from the fragmentation of a color-neutral system of par- 
tons. In these systems correlations are produced by local 
conservation of charge, flavor, energy and momentum in 
the strong interaction, and by quantum statistics. In 
high energy heavy ion collisions aspects of these elemen- 
tal correlations might persist, especially at high trans- 
verse momentum (pt) since the "memory" of the early 
hard partonic scattering is not easily erased there. In 
contrast, minijets at lower pt are expected to have 
shorter mean free paths in the medium and thus are 
more likely to dissipate, erasing correlations. The col- 
lision overlap density and size of the interaction volume 
are changed by varying the centrality, which might also 
control the degree of equilibration in these systems. We 
study the correlation structure in peripheral collisions, 
caused by minijets, which evolves with centrality and pt 
in a manner suggesting strong dissipation of minijet frag- 
mentation by the longitudinally expanding medium. 

The data presented here were obtained with the STAR 
Time Projection Chamber (TPC) [3, mounted inside a 
solenoidal magnet. Charged-particle tracking with the 
TPC covers large acceptance well suited for precision 
studies of correlation structures over a wide range of 
scales. The minimum-bias event trigger discriminates on 
a neutral-spectator signal in the Zero Degree Calorime- 
ters 0. Central events were selected by additionally re- 
quiring a high charged particle multiplicity within |?7| < 1 
in the Central Trigger Barrel scintillators 0- Accepted 
charged-particle tracks had > 15 TPC space points and 
>52% of the estimated maximum possible number of 
space points (to eliminate split tracks), passed within 
3 cm of the event vertex and were within the kine- 
matic acceptance: \ri\ < 1, full 2tt in azimuth, and 
0.14 < pt < 2.1 GeV/c. Accepted events had their 
primary vertex within 25 cm of the geometric center 
of the TPC longitudinally and had >15 accepted TPC 
tracks. About 0.6 M central and 0.3 M peripheral events. 



recorded in the y^s^vw = 200 GeV run, were analyzed. 

Two-point correlations and power spectra of point-to- 
point fluctuations are complementary measures used to 
study the correlation structure of random fields (such 
as TPC events). The former has computational com- 
plexity 0{N'^) {N is event multiplicity). The latter, 
implemented via the discrete wavelet transform (DWT) 
method, is 0{N). The DWT-based dynamic texture 
measure, defined below, is used in this work and was 
originally applied to relativistic Pb-|-Pb collisions by 
NA44 113. 

In this approach, the measured particle distribution 
p{(f>, rj) in a single event is expanded in the complete 
orthonormal wavelet basis of Haar The scale of 

this basis is defined by the scaling function g{x) ~ 1 
for < X < 1 and otherwise. The function f{x) — 
{1 for < a; < 0.5; -1 for 0.5 < x < 1; else 0} is the 
wavelet function. The experimental acceptance in 77, (/), 
and Pt is split into equal bins in 7^,0 and pt bins ex- 
ponentially growing to equalize bin statistics. To keep 
notation simple but explicit, we introduce rj' = (77 -I- l)/2 
and (j)' = 0/27r so that 77', </>' £ [0, 1]. The scaling function 
of the Haar basis in two dimensions G{4>, 77) = g{4>')gi'>l') 
is just the bin acceptance (modulo units). The wavelet 
functions F"*" (the directional sensitivity mode A is ei- 
ther along azimuth (j), pseudo-rapidity 77, or diagonal 
^7; directions) are F'f"^ = f{(j}')f{i), F't' = f{(l)')g{i), 
F"^ = g{4'')f{'ri'). We define a two dimensional wavelet 
basis: 

F,^„^,_^.(</),7?) =2-F^(2V-*,2'"77'- j), (1) 

where m > is the integer scale fineness index [T^ , inte- 
gers i and j index the positions of bin centers in (p' and 
77', and < i,j < 2"*. Scaling functions Gm,i,j{(t>iV) are 
constructed analogous to EqQ] Arbitrary density p((/), 77) 
is expanded as 

p((/),77) = (p,Go,o,o)Go,o,0+ iP^Frn,^J)F^,^,J, (2) 

m,i,j.X 

where (p, G) and (p, F^) are expansion coefficients ob- 
tained by projecting density p(0, 77) onto the basis func- 
tions. 

In practice m < TTimax, where rrimax is the finest scale 
limited by track resolution and, due to the needs of event 
mixing, by the number of available events. The coarser 
scales correspond to successively re-binning the track dis- 
tribution. The analysis is best visualized by considering 
the scaling function Gm.i,ji4'i v) as binning the track dis- 
tribution p(0, rj) in bins i,j of given fineness m, while the 
wavelet expansion coefficients {p,Fj^ ^ j) give the differ- 
ence distribution for data with binning one step finer. 
The wavelet expansion coefficients were calculated using 
the code waili [T3l |. 

The power spectrum is defined as 



P\m)^2-'"^Y{p,F^l,^^y, (3) 
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where the overhne denotes an average over events. 
P^irn) is independent of m for an uncorrelated p. How- 
ever, for physical events depends on m due to the 
presence of static texture features such as acceptance 
asymmetries and imperfections (albeit minor in STAR), 
and non-uniformity of dN / drj. To remove these known 
features from the analysis a reference is constructed from 
mixed events starting with individual (i^, ?/) pixels of true 
events at the finest scale used in the analysis (16 x 16). A 
"mixed event" consists of 16 x 16 {(j), rj) pixels from true 
events, where each pixel is taken from different, but sim- 
ilar, real events. The power spectrum P^^^ is obtained 
from Eq. lousing the expansion coefhcients in Eq.|21 f'^ix 
contains static, experimental track density artifacts plus 
statistical noise. The quantity of interest is the differ- 
ence, P^y^ = -Ptruc ~ ^mix' Called dynamic texture [lo| . 

In studying the dynamic texture data as a function of 
Pt, the desirable normalization is such that the results 
are independent of pt bin size under the assumption of 
large-scale correlations in pt {i.e. larger than the pt ac- 
ceptance). In this case for increasing number of particles 
N in an increasing pt bin, i^dyn °^ while P^i^ oc iV, 
being a Poissonian variance. Therefore we present the 
data as the combined quantity -PdynZ-^mix/-^- 

Systematic error can be introduced in P^y^ by the pro- 
cess of event mixing. For example, events with different 
vertex positions along the beam axis are reconstructed 
with slightly different efhciencies and acceptances with 
respect to ry. This variable efficiency may fake a dynamic 
texture effect in rj. In order to minimize such errors, 
events were grouped into event classes with similar mul- 
tiplicity (within 50) and vertex position (within 10 cm), 
^dyn constructed using only events within each of 
these two classes. Results showed no vertex dependence. 
The upper limit on the systematic error due to z-vertex 
position variation is set by the statistical error of the 
data, shown in the figures. 

Event centrality in this analysis is characterized by the 
accepted number of quality tracks in the TPC and ex- 
pressed a s a p ercentage of the total inelastic cross-section, 
as beforeflj. Event classes in multiplicity are grouped 
to form two centrality classes: central^ with 4% of the 
most central events, and peripheral, with event centrality 
varying between 60% and 84%. The hijing gener- 
ator events for the Monte Carlo comparison are selected 
to match these centrality ranges. 

Track splitting (one particle reconstructed as > 1 
track) contributions were eliminated by track quality re- 
quirements. Track merging (> 1 particles reconstructed 
as one track) mocks up anticorrelations and can induce 
systematic error. To estimate this effect, central HI- 
JING events were filtered with an algorithm emulating 
track recognition properties of the TPC (1^]. The sim- 
ulation results can be expressed as a set of coefficients 
relating Ptruc, Pmix and dN/ dpt in the original and fil- 
tered HIJING data. An estimate of track merging effects 
in the data was obtained from the inverse of these co- 
efficients. The resulting systematic error was estimated 
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FIG. 1: Peripheral events; normalized dynamic texture for 
fineness scales m = 0,1,0 from left to right panels, respec- 
tively, as a function of pt. • - STAR data; solid line - hijing 
without jet quenching; dash-dotted line - hijing without jets. 



to be 0.5 X 10^^. Systematic error due to non-primary 
background was estimated assuming that the correla- 
tions between true primary and non-primary particles 
could be anything from zero to that of primary particles 
themselves. The systematic error estimate was taken to 
be half the difference between these two limits which is 
10% of the signal at pt ^ 0.2 GcV/c, falhng to 3.5% at 
Pt = 1 GeV/c. This estimate applies to both centrality 
classes. 
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FIG. 2: Central events: normalized dynamic texture for fine- 
ness scales m = 0, 1,0 from left to right panels, respectively, 
as a function of pt. • - STAR data; solid line - hijing with- 
out jet quenching; dashed line - hijing with quenching; □ - 
peripheral STAR data from Fig. prenormalized as described 
in the text. The rectangles around two chosen points show 
the estimated systematic errors. 

Fig. n presents measured large-scale dynamic texture 
in peripheral collisions compared with hijing predictions 
where fineness scale m = is used for the rjcji and 77 
modes and fineness to = 1 is used for the (j) mode. The 
finer scale for Scf) is used so that the angular coverages 
subtended by 6(j) for to = 1 and Srj for to = are com- 
parable. This scale, with (j)-hm size S(j) = 7r/2 is 
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dominated by elliptic flow The HIJING calculations 
without jet quenching show a region of approximately 
constant signal near pt ~ 0.5 GeV/c followed by an in- 
crease for Pt > 0.8 GeV/c, obtained by "turning on" jets 
in the model. In that pt range the STAR data also in- 
crease with Pt- Momentum conservation suppresses the 
difference in the numbers of tracks emitted in the oppo- 
site directions. This effect is absent in the mixed events, 
resulting in negative P^y^, seen in when jets in hijing 
are "off" . Comparing the two simulations in Fig. ^ we 
conclude that fluctuations in local hadron density due 
to jets are observable in peripheral RHIC collisions at 
0.8 < Pt < 2 GeV/c. This supports but does not prove 
the identification of similar signals in the data at these 
Pt with minijets. Without ruling out other sources of 
angular correlations at such pt, we use Occam's razor 
to adopt the well established effect - fragmentation of 
semi-hard scattering products (jets or minijets) - as the 
explanation. 

Central event data and hijing predictions with and 
without jet quenching arc shown in Fig. |21 The most 
striking difference here compared to the peripheral data 
in Fig. n is the reduction in the magnitude of the P^y^ at 
larger pt > 0.6 GeV/c, the data becoming slightly neg- 
ative near 1 GeV/c in sharp contrast to the jet-like be- 
havior predicted by hijing. The perturbative partonic 
energy loss model of jet quenching in hijing seems to 
miss the correlation aspect of the picture, at least at 
these Pt- In the absence of a successful theory to de- 
scribe the effect, we formulate and test a "null hypothe- 
sis": the correlation structure ^'dyn/^mix Au+Au col- 
lisions is independent of centrality. Then, the difference 
in ^dynZ-^mix/-^ Central and peripheral events (includ- 
ing the Pt trends) is due to the difference in 1/N (i.e. in 
dN/dpt, for N = Nipt) = dNipt)) M- Shown in 

Fig. 121 as symbol □ are the peripheral data from Fig. ^ 
rescaled under an assumption of the "null hypothesis" by 
xA^(pi)|po,iph/A^(pt)|conti- The left panel shows that the 
7;0-mode is less affected by centrality, reflecting a super- 
position of the opposite centrality trends in ij and ip. We 
hypothesize that the deviation of the STAR data from 
the "null hypothesis" in t] in the otherwise correlated 
system points to a randomization (dissipation) of mini- 
jet structure in the longitudinal direction. Longitudinal 
expansion of the hot, dense medium formed early in the 
collision singles out the rj direction and is likely to be part 
of the dissipation mechanism. If so, at pt > 0.6 GeV/c we 
may be observing an effect of the longitudinally expand- 
ing medium on parton fragmentation or hadronization. 

In each panel of Figs.^andl^lthe dynamic texture data 
increase with decreasing pt for pt < 0.4 GcV/c. Data 
stay non-zero at low pt for all three modes in the experi- 
ment and for the ?7-mode in hijing. In this pt range, the 
correlations are likely dominated by centrality-dependent 
effects such as the final state quantum statistical inten- 
sity interference. Coulomb effect and longitudinal string 
fragmentation physics, simulated in hijing. Modifica- 
tion of the latter effect with centrality is the subject of a 
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FIG. 3: Scale dependence of the dynamic texture in peripheral 
and central events. (a,c): 0.2 < pt < 0.28, (b,d): 1.1 < pt < 
1.5 GeV/c. • - STAR data; solid line - hijing without jet 
quenching; dashed line - hijing with quenching. A systematic 
error estimate is shown as a hatched area. Errors on different 
scales are estimated independently. 



separate publication[l9j. 

Fig. 131 shows a scale dependence of the ry-mode in the 
low and higher pt intervals. At low pt, the peripheral 
and central trends qualitatively agree, whereas at higher 
Pt, a modification with centrality is seen, which testi- 
fies to the presence of new physics at higher pt. The 
reduction of the dynamic texture in central events with 
respect to both HUING and the peripheral STAR data 
is most dramatic at the coarser scales. The longitudinal 
expansion correlates 77 with the longitudinal coordinate 
z, and z - with time. Final state particles with large 
6r] arc more likely to be separated by a space-like inter- 
val. Thus, the larger Sr] correlations are more likely to 
have their cause in the particles' common past, reflect- 
ing the early stage of the system, whereas the fine scale 
features are formed later under conditions little differ- 
ent from peripheral collisions or conventional hadronic 
models. The negative P^y^^ in Fig. EJd) points to the 
presence of an anticorrelation mechanism, which could 
include existence of a characteristic scale in the longitu- 
dinal separation of hadrons in the course of hadroniza- 
tion. Lack of scale dependence in Fig. |2Id), relative to 
Fig. |3Jb), may be contrasted with progressive reduction 
of small-scale Fourier harmonics from hadronic diffusion 
discussed in |^ . Alternatively, pre-hadronic transport on 
77 involving partonic diffusion could provide a more effi- 
cient equilibration mechanism. Other mechanisms such 
as convective turbulent transport |20l | might also play a 
role. The reduction of dynamic texture reported in this 
Letter provides a new quantitative argument in favor of 
equilibration or dissipation effects. However, we observe 
that the hadronic final state is not correlation free, even 
for central events. 

In summary, a non-trivial picture emerges when the 
DWT power spectrum technique is applied for the first 
time to Au4-Au collision data from RHIC. Large-scale 
[St] = 1) angular correlations for pt < 2.1 GeV/c are 
observed in peripheral events and identified with mini- 
jets. In central events, those correlations are suppressed 
with increasing pt and Sr]. This indicates a major change 
in the properties of the medium with increasing colli- 
sion centrality, implying the development of a dissipative 
medium. In the course of its longitudinal expansion, this 
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hypothetic medium influences via interactions the struc- 
ture of correlations, inherited from the kinematics of the 
initial-state semi-hard scattering, causing their dissipa- 
tion and partial equilibration. 
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